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ABSTRACT 

Among the many factors that determine overall rocket performance, propellant density is important because it affects the size 
of the rocket. Thus, in order to decrease the size of a rocket, it may be desirable to increase the density of propellants. This 
study analyzes the concept of increasing the propellant density by employing a cooling source submerged in the liquid 
propellant. A simple, mathematical model was developed to predict the rate of densification and the propellant temperature 
profile. The mathematical model is generic and applicable to multiple propellants. The densification rate was determined 
experimentally by submerging a cooling source in liquid oxygen at constant, positive pressure, and measuring the time rate of 
change in temperature with respect to vertical position. The results from the mathematical model provided a reasonable fit 
when compared to experimental results. 


1. Introduction 

In an attempt to lower the inert mass of a launch 
vehicle, propellant densification is viewed as a method to 
reduce the propellant tank volume and a means to reduce 
the mass of the propellant tank and, therefore, increase 
payload mass. Many have conducted studies for propellant 
densification [1] [2] [3]. Most of the studies suggest to 
place the densification operation between the ground 
storage tank and launch vehicle [4] [5] [6]. This paper 
investigates placement of the densification operation 
within the ground storage tank. Also, this paper predicts 
the rate of propellant densification using a cooling source 
within the liquid propellant and compares the prediction 
to experimental results. This work is similar to the 
Integrated Refrigeration and Storage method proposed by 
Notardonato [7] [8]. 

2. Material and methods 
2.1 Test Article 

The main test article is comprised of a 400L dewar 
with a 0.254 m neck opening that accepts a male bayonet. 
Eden Cryogenics constructed the dewar using 304 
stainless steel with an inside diameter of 0.76 m and a 
maximum allowable working pressure of 515 kPa 
absolute. The male bayonet, also constructed with 304 
stainless steel, incorporated two fluid penetrations and 
two instrumentation cable penetration. A heat exchanger 
extended below the male bayonet and connected to the 
male bayonet by 0.013m and 0.025m stainless steel tubing, 
which insulated with a thickness of 0.076m of 
polytetrafluoroethylene. The fluid penetrations allowed 
coolant to enter the heat exchanger through the 0.013m 
diameter tubing and exit the heat exchanger through the 
0.025m tubing. The heat exchanger was 0.013m diameter 
copper tubing with annular fins totaling 0.36m in length. 

Peripheral equipment included a liquid nitrogen 
source, liquid nitrogen subcooler, and gaseous oxygen 
source. Liquid nitrogen, in 110L and 180L dewars 
saturated at 0.21 kPag, provided the cooling fluid for the 


heat exchanger inside the main test article. The liquid 
nitrogen flowed through a simple liquid nitrogen 
subcooler to ensure the liquid nitrogen is subcooled before 
entering the main test article. The liquid nitrogen 
subcooler consisted of copper heat exchanger submerged 
in liquid nitrogen subjected to atmospheric pressure. The 
gaseous oxygen source, from a 43L k-bottle at ambient 
temperature and an initial pressure of 16.6 kPag, provided 
the means to control the pressure within the main test 
article. 

2.2 Instrumentation 

A micarta instrumentation rake was attached to the 
bottom of the neck flange. A total of 28 RTD measured the 
temperature inside the test dewar. Twenty RTD measured 
the temperature of position in vertical and the remaining 
eight RTD measured the horizontal temperature at two 
vertical positions. One pressure transducer measured the 
test dewar pressure, while a separate pressure transducer 
measured the pressure in the coolant system. Two 
thermocouples measured the temperatures of the coolant 
entering and exiting the male bayonet. Two mass flow 
controllers measured and controlled the quantity of 
coolant flowing through the heat exchanger and the 
quantity of gaseous oxygen flowing into the test dewar. 
The specifications for the nitrogen and oxygen were MIL- 
PRF-27401F Grade B and MIL-PRF-25508G Grade F, 
respectively. 
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Figure 1: Mechanical and Instrument Schematic of Test 
Materials 


2.3 Operation 

The test dewar was filled with 200L of liquid oxygen, 
so that the liquid height above and below the heat 
exchanger is 0.152m and 0.305m, respectively. Liquid 
nitrogen is used as the coolant and flowed from the liquid 
nitrogen dewar through the heat exchanger inside the test 
dewar and to the nitrogen mass flow controller. The 
nitrogen mass flow controller was set at its maximum 
value of 100 standard liters per minute. The gaseous 
oxygen source maintained constant pressure within the 
test dewar, by flowing through the oxygen mass flow 
controller. The gaseous oxygen mass flow controller was 
controlled using a PID controller within Labview 
referencing the main test article pressure, P-107. Three 
tests were run at a constant pressure of 3 psig, 5 psig, and 
7 psig. 

3. Mathematical Model 

Prior to activating the internal heat exchanger, the 
steady-state condition of the main test dewar begins with a 
thermally stratified ullage with respect to vertical position 
and a uniform temperature throughout the liquid. The 
temperature of the liquid is saturated temperature at the 
main test dewar storage pressure. Once the internal heat 
exchanger is activated, the liquid is divided into two zones 
depending on whether the liquid is above or below the 
internal heat exchanger. A mathematical model is 
provided for both zones. 



Figure 2: Main Test Article at Steady- State, Initial 
Condition 



Figure 3: Main Test Article During Heat Exchanger 
Operation 

3.1 Liquid Zone Above Heat Exchanger 

The liquid zone above the heat exchanger is modeled 
using the heat equation with constant temperature 
boundary conditions at both, the top condition and bottom 
condition. The top boundary is the vapor-liquid interface, 
which is always at the saturated temperature at the system 
pressure. Because the main test dewar is operated at 
constant pressure, the temperature at the top boundary 
remains constant. The temperature at the bottom 
boundary is given by the temperature of the internal heat 
exchanger. The heat equation (1) is solved using the 
boundary and initial conditions (2), (3], and (4] to arrive at 
the modeled temperature change rate for the liquid zone 
above the heat exchanger (5). 
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3.2 Liquid Zone Below Heat Exchanger 

The liquid zone below the heat exchanger is modeled 
using the heat equation averaged across the entire zone 
height to estimate the convective heat transfer. The upper 
boundary condition, modeled as a constant temperature 
boundary, is given by the temperature of the internal heat 
exchanger. The lower boundary condition is taken to be 
insulated. The liquid is initially uniform in temperature 
with respect to vertical height. The heat equation (6) is 
solved using the boundary and initial conditions (7), (8], 
and (9). The solution is provided by (10). 


Governing Equation: 
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Because the temperature of the heat exchanger is always 
cooler than or equal to the temperature of the liquid 
oxygen below the heat exchanger, natural convection 
occurs in the liquid zone beneath the heat exchanger. The 
natural convection provides a uniform temperature with 
respect to height at all times. Thus, the temperature 
profile, (10), is averaged over the entire height of the 
liquid zone beneath the heat exchanger to arrive at the 
modeled temperature change (11). 
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4. Results 

4.1 Test Results for Liquid Zone Above Heat Exchanger 
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Figure 4: Test Results for Liquid Zone Above the Heat 
Exchanger (3 psig) 
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Figure 5: Test Results for Liquid Zone Above the Heat 
Exchanger (5psig) 
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Figure 6: Test Results for Liquid Zone Above the Heat 
Exchanger (7psig) 

4.2 Test Results for Liquid Zone Below Heat Exchanger 
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Figure 7: Test Results for Liquid Zone Below the Heat 
Exchanger (3psig) 
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Figure 8: Test Results for Liquid Zone Below the Heat 
Exchanger (5psig) 
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Figure 9: Test Results for Liquid Zone Below the Heat 
Exchanger (7psig) 


5. Discussion 


The results from Figure 4, 
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Figure 5, and Figure 6 show that the liquid zone above 
the heat exchanger experiences conductive heat transfer. 
Minimal convection occurs within the liquid zone above 
the heat exchanger as shown by the stratification that 
grows throughout the heat exchanger operation. The 
mathematical model, Equation (5), predicts the 
temperature change rate at the T5 position with a slight 
lag. However, with exception of the results in Figure 4, the 
lag increases for the T6 position prediction in temperature 
change rate. 

The results from Figure 7, Figure 8, and Figure 9 show 
the liquid zone below the heat exchanger experiences 
convective heat transfer. All four position temperature 
sensors, T1 through T4, maintain the same temperature 
throughout the heat exchanger operation, which implies a 
fluid well mixed by convection currents. The mathematical 
model, Equation (11), predicts a temperature change that 
is less than the experimental results. The reason for the 
temperature change rate lag could be that Equation (11), is 
an average for a conduction model and does not model the 
convection currents, which increases the heat transfer 
from the upper boundary at the internal heat exchanger. 

The mathematical models presented neglect the heat 
flux that enters the dewar from the sides in Equation (5), 
and the heat flux that enters the dewar from the sides and 
bottom in Equation (11). Thus, Equation (5) and Equation 
(11) produce models that are one-dimensional. However, 
from the horizontal temperature sensors, T21 through 
T24, showed no gradient in the radial direction and was 
the same temperature as the centerline temperature at the 
same vertical height. Therefore, the one-dimensional 
mathematical models are reasonable. 


Because the mathematical model for the liquid zone 
beneath the heat exchanger is an average temperature 
throughout the liquid zone, the temperature at the upper 
boundary does not equal temperature of the heat 
exchanger. Thus, a temperature discontinuity exists 
between the two models, Equation (5) and Equation (11), 
at the heat exchanger, and the possibility exists for colder 
fluid positioned above warmer fluid. In reality, natural 
convection occurs when a colder, more dense fluid is 
above a warmer, less dense fluid. Therefore, the 


temperature of the liquid above the heat exchanger can 
never be less than the temperature below the heat 
exchanger. Thus, the discontinuity between (5] and (11] 
can be corrected by mandating that temperature of the 
liquid above the heat exchanger be the maximum of either 
Equation (5] or Equation (11). Equation (12] and 
Equation (13] provide the unified temperature model. 
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T(y,t) = Equation (11) for 0 < y < hi (13) 


Figure 10 shows the predicted vertical temperature 
profile, during heat exchanger operation. The vertical 
temperatures are represented by T1-T12, which are 
located vertically in the main test dewar. The vapor-liquid 
interface is at the T8 and the heat exchanger is located at 
the T4. The left portion of Figure 10 shows the 
temperature discontinuity between the mathematical 
models, (5) and (11). The right portion of Figure 10 shows 
the unified temperature models from Equation (12) and 
Equation (13). Figure 11 shows the combined 
experimental data of Figure 6 and Figure 9 for comparison 
with Figure 10. 
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Figure 10: Exposure and Correction of the Discontinuity 
Between Both Mathematical Models 
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Figure 11: Test Results for the Liquid (7psig) Shown as 
the Temperature Profile for the Vertical 
Position. 

The coolant removed about 330W from the main test 
article. However, not all of the heat removed was used to 
cool the liquid oxygen. The heat transfer from the ambient 
was approximately 20 W, according to the experimental 
normal evaporation rate. Also, a portion of the gaseous 
oxygen that flows into the ullage is cooled and liquefied, 
accounting for 180W. Therefore, 130W of heat was 
removed from the liquid oxygen. 

The major benefit for densification within the storage 
tank under positive pressure allows the use of a low power 
cooling source that densities the propellant over time 
while inhibiting atmospheric intrusion. Same species 
pressurization also assures that propellant purity will be 
maintained throughout the densification process. 
However, as shown above, the liquid below the heat 
exchanger can be fully densified, while the liquid above the 
heat exchanger is partially densified. Also, cooling the 
pressurant gas required about 55% of the cooling available 
from the heat exchanger, decreasing the densification rate 
of the liquid. 

6. Conclusions 

The use of a submerged cooling source is a viable method 
to obtain densified and sub-cooled propellants. Also, same 
species pressurant can be used to maintain positive 
pressure of a storage tank during the densification 
operation. This study confirmed that conduction heat 
transfer dominates the fluid above the heat exchanger, 
while convection heat transfer dominates the fluid below 
the heat exchanger. 
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8. Variables 

n integer used for summation ( ) 

h height (m) 

t time (s) 

T temperature (K) 

r non-homogeneous component of heat equation 

y vertical component of coordinate system (m) 

a thermal diffusivity (m * 1 2 /s) 

Subscripts 

1 vertical height at heat exchanger 

2 vertical height at vapor liquid interface 

sat saturated temperature 

s heat exchanger surface temperature 


